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Abstract 



A search has been made for massive resonances decaying into a quark and a vector 
boson, qW or qZ, or a pair of vector bosons, WW, WZ, or ZZ, where each vector boson 
decays to hadronic final states. This search is based on a data sample corresponding to 
an integrated luminosity of 5.0 fb^^ of proton-proton collisions collected in the CMS 
experiment at the LHC in 2011 at a center-of-mass energy of 7 TeV. For sufficiently 
heavy resonances the decay products of each vector boson are merged into a single 
jet, and the event effectively has a dijet topology. The background from QCD dijet 
events is reduced using recently developed techniques that resolve jet substructure. 
A 95% CL lower limit is set on the mass of excited quark resonances decaying into 
qW (qZ) at 2.38 TeV (2.15 TeV) and upper limits on the cross section for resonances 
decaying to qW, qZ, WW, WZ, or ZZ final states. 



Submitted to Physics Letters B 



*See Appendix A for the list of collaboration members 



1 



1 Introduction 

New resonances that decay preferentially into hadronic final states are of particular interest in 
a variety of scenarios for physics beyond the standard model (SM) HHH. Searches for events 
having a pair of hadronic jets with large invariant mass have been performed by the Compact 
Muon Solenoid (CMS) and ATLAS experiments at the Large Hadron Collider (LHC) [10, 11 J. In 
principle, these searches are also sensitive to final states that include one or two massive vector 
bosons W/Z because the vector bosons have large hadronic branching fractions and because 
their masses are much smaller than those of the hypothetical parent states, i.e. they are highly 
"boosted". This implies that the pairs of quark jets produced by the vector boson decays merge 
into single W/ Z-jets in a real detector. The sensitivity of the present large mass dijet searches is 
limited by the presence of background from ordinary strong interaction processes that produce 
pairs of quark and gluon jets. 

The analysis presented here exploits the enhancement of the sensitivity of a standard dijet anal- 
ysis for processes that produce W/Z-jets in the final state by the application of techniques that 
can identify W/Z-jets and suppress quark and gluon jets ("W/Z-tagging"). This CMS study 
is performed on pp collision data at a center-of-mass energy of 7 TeV, corresponding to an in- 
tegrated luminosity of 5.0 fb^^. We consider events with two high-transverse-momentum jets 
in the final state. We identify "subjets" inside jets using recent developments in the area of jet 
substructure [12J. Pairs of subjets are used to explicitly reconstruct W or Z bosons, therefore 
substantially suppressing backgrounds from quantum chromodynamics (QCD) interactions. 
This search follows closely the procedures of the corresponding dijet search [10], performed in 
the same dataset, but with strongly reduced QCD background because of the W/Z-tagging. 

We consider three benchmark scenarios that would produce singly- or doubly-tagged events: 
an excited quark q* [4J decaying into a quark and a W or Z boson; a Randall-Sundrum (RS) 
graviton decaying to WW or ZZ Grs [Bl; and a heavy partner of the SM W boson W which 
decays to WZ fSI. The most stringent limits on the q* model have been set in dijet resonance 
searches at the LHC by considering the qg final state fW] or inclusively all-hadronic final 
states [11]. The most stringent lower limit (at 95% CL) on the q* mass to date is 3.3 TeV [lOJ. 
Specific searches for the qW and qZ final states have previously been reported at the Teva- 
tron |[14[ |T5|, which exclude resonances decaying to qW or qZ with masses up to 540 GeV, 
and at the LHC fl^, which extends the mass exclusion of qZ resonances up to 1.94 TeV. For 
the Grs, there are phenomenological models favoring the decay of the Grs into vector bosons 
rather than photons or fermions |[T7] - [T9l . In particular, the ZZ final state has been explored 
experimentally l|20H22) . setting lower limits on the Grs mass as a function of the coupling pa- 
rameter k/ Mpi, where k is the curvature of the warped space and Mpi the reduced Planck mass 
(Mpi = Mpi/ y/Sn). For the W , the most stringent limits are reported in searches with leptonic 
final states ||23[|24| , and the current lower limit on the W mass is 2.5 TeV. The limit varies by 
0.1 TeV, depending on the chirality of the W couplings. Specific searches in the WZ final state 
have also been reported [.25n26il setting a lower limit of 1.1 TeV. 

This paper is organized as follows. First, the CMS detector, and the simulated and collision data 
samples on which the analysis is based are briefly described. Then the event reconstruction, 
and selection are detailed, and the W/ Z-tagging technique is described. The following section 
describes the modeling of detector acceptances and signal efficiencies as well as the validation 
of the W/Z-tagging techniques with data. After this follows the description of the modeling 
of the background, the systematic uncertainties and the limit setting procedure. Finally the 
results and conclusions are presented. 
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4 Event reconstruction and selection 



2 CMS detector 

The CMS detector II27I is well-suited to the reconstruction of hadronic jets because it incor- 
porates finely segmented electromagnetic and hadronic calorimeters, and a charged-particle 
tracking system. Charged particles are reconstructed in the inner tracker, which is immersed 
in a 3.8 T axial magnetic field. The inner tracker consists of three cylindrical layers and two 
endcap disks at each end of silicon pixel detectors, and ten barrel layers and twelve endcap 
disks at each end of silicon strip detectors. This arrangement results in full azimuthal coverage 
{0 < (p < 2n) within |^| < 2.5, where i] is the pseudorapidity defined as t] = — ln[tan(0/2)]. 
CMS uses a polar coordinate system, with the z axis coinciding with the beam axis; 6 is the 
polar angle defined with respect to the positive z axis. Muons are measured in gas-ionizing de- 
tectors embedded in the steel return yoke. A lead-tungstate crystal electromagnetic calorimeter 
(ECAL) up to |f/| = 3 and a brass /scintillator hadronic calorimeter (HCAL) up to |f/| =5 sur- 
round the tracking volume and allow photon, electron, and jet reconstruction. The ECAL and 
HCAL cells are grouped into towers projecting radially outward from the interaction region. 
In the central region {\rj\ < 1.74) the towers have dimensions At] = A(p = 0.087; at higher rj, 
the Ai] and widths increase. For optimum jet reconstruction, the tracking and calorimeter 
information is combined in an algorithm called particle flow [28J, which is described below. 

3 Simulated and collision data samples 

The sample of proton-proton collision data at i/s = 7TeV, corresponding to an integrated 
luminosity of 5.0 fb^^, was collected in 2011. The events were collected using the logical "or" 
of a set of triggers based on requirements on Hj = ^gjg Ej (Ej is the transverse energy of a 
jet) and the invariant mass of the two highest pj jets in an event, whose thresholds were raised 
progressively to cope with an increase in the peak luminosity during 2011. 

Data are compared to Monte Carlo (MC) simulations of the QCD background generated using 
both PYTHIA6.424 [29J and HERWIG++ 2.4.2 [30]. PYTHIA 6 is used with CTEQ61L |31J and 
HERWIG++ with MRST2001 [32J parton distribution functions. Tune Z2 ||33l iMl is used with 
PYTHIA 6, while the tune version 23 EOl is used with HERWIG++. 

The sensitivity of the event selection to the benchmark processes is evaluated using simulated 
samples of events from excited quarks, RS gravitons, and W production and decay models. 
The process q* — > W/Z + jet is generated using PYTHIA 6. The process Grs — > WW/ZZ is 
generated using HERWIG++ and its cross section is taken from PYTHIA 6. While HERWIG++ 
contains a more detailed description of the angular distributions than PYTHIA 6 for this pro- 
cess [1351 ■ the cross section is taken from PYTHIA 6 which has been used as a reference model in 
related analyses fTO]. RS graviton production is studied with k/ Mpi = 0.1, which determines a 
resonance width of about 1% of the resonance mass which is about a factor 5 smaller than the 
experimental resolution for dijets. This width is much smaller than suggested by the model in 
Ref. fl 7| , which predicts resonance widths of the order of the experimental resolution, allowing 
for interpretation in this model only approximately. The process W — )■ WZ is generated using 
PYTHIA 6 with Standard Model V — A couplings and without applying k-factors. All Monte 
Carlo events are passed through the CMS detector simulation based on Geant4 II36II . 

4 Event reconstruction and selection 

Events are reconstructed using the particle flow algorithm, which attempts to identify and 
measure all the stable particles in a collision by combining information from all the subdetec- 
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tors. This algorithm categorizes all particles into five types: muons, electrons, photons, charged 
and neutral hadrons. The resulting particle flow candidates are passed to the anti-A;T l|37l and 
Cambridge-Aachen (CA) ||38ll39| jet clustering algorithms, as implemented in FastJet fiO'.ITl 
to create jets. A distance parameter of size R = ^(Af/)2 + (A(/))2 = 0.8 is used for the CA 
algorithm, while R = 0.5 is used for the anti-Zcx algorithm. While the anti-A;x jets are used 
to select events and reconstruct the dijet invariant mass mjj, the CA jets are used to identify 
events containing hadronically-decaying W or Z bosons. This choice has been made because 
the CA algorithm was found to be more efficient (for the same mistag rate) at finding hard sub- 
jets within the jets in simulation-based studies ||T2|, while the anti-fcx jets have the best energy 
calibration. 

Events must have at least one reconstructed vertex within |z| < 15cm, to suppress back- 
grounds solely triggered by calorimeter noise. The primary vertex is defined as the vertex 
with highest sum of squared track transverse momenta (pj). Charged particles not originating 
from it are removed from the inputs to the jet clustering algorithms. This requirement removes 
particles which arise from additional pp interactions in the same pp bunch crossing (pileup 
interactions). An event-by-event jet-area-based correction |'42'-'44j is applied to remove the re- 
maining pileup energy which is due to neutral particles originating from the other vertices. The 
pileup-subtracted jet four momenta are finally corrected to account for the difference between 
the measured and true responses to hadrons [44J. When jets are decomposed into subjets, as de- 
scribed later, the energy estimate relies on the calibrated reconstructed input particles without 
further corrections. 

Events are initially selected by requiring that they have at least two ardi-kj jets with px > 
30 GeV and \ < 2.5. The two highest-px jets are required to have a pseudorapidity separation 
|A7/| < 1.3 to reduce the QCD dijet background l|45| . Finally, the dijet invariant mass is required 
to be larger than 890 GeV. This threshold is defined by the triggers, which were found to be 
99% efficient for dijet events with masses above this threshold. 

In events passing this selection, "boosted" (high px) hadronically-decaying W or Z bosons are 
identified with a W/Z-tagging algorithm using jet pruning [46J, a technique which removes 
the softest components of the jets. In the jet pruning technique HZlllHl, a jet is reclustered using 
all the particles used to build a CA jet, ignoring in each recombination step the softer "protojet" 
if the recombination is softer than a given threshold z^ut or forms an angle AR wider than Dcut 
with respect to the previous recombination step. The hardness of a recombination z is defined 
as z = mm{pj, pj)/ pj, where pj and pj are the px of the two protojets to be combined and pj is 
the Px of the combined jet. The following selection is then applied to the pruned jets to identify 
jets from hadronic W/Z decays by exploiting the variables used in Ref. [49 1. The total pruned 
jet mass mjet must satisfy 70 GeV < mjet < 100 GeV. Two subjets are obtained by reversing the 
last clustering iteration of the pruned jet clustering. The mass drop ^ is defined as the ratio of 

the heaviest subjet mass (defined as mi) and the pruned jet mass. The mass drop is required to 
satisfy ^ < 0.25, which eliminates the soft divergent components in QCD jets. These criteria 

are designed to select W and Z candidates in which the subjets are similar in energy and mass. 

Comparisons of the dijet invariant mass distributions for untagged, single-tagged, and double- 
tagged event samples are shown in Fig. [l] The data are shown as solid points and the PYTHIA 
6 and HERWIG++ simulations are shown as solid red and dashed blue curves, respectively. 
The simulations are normalized to the number of data events in each category and the shapes 
are compared; the agreement of the normalization driven by the W/Z-tagging efficiency is 
discussed in the next section. The pythia 6 prediction is found to agree with the data while the 
HERWIG++ prediction decreases more steeply with mass. However, no systematic uncertainties 
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5 Signal characterization 
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Figure 1: Comparisons of the dijet invariant mass distributions between data and Monte Carlo 
(PYTHIA 6 and HERWIG++) simulations. The three sets of lines correspond to the inclusive dijet 
category (no W/Z-tag required), single W/Z-tagged, and double W/Z-tagged events. The 
simulations are normalized to the number of data events in each category. 

are taken into account and only the dominant background from QCD interactions is considered. 
In this analysis, the background shape is modeled from the data themselves. Therefore, the 
analysis depends on QCD simulation only to provide guidance and cross checks. 



5 Signal characterization 

A search for dijet resonances corresponding to several benchmark physics models is performed. 
Using the W/Z-tagging algorithm, both single W/Z-tag and double W/Z-tag events are ex- 
amined. The signals that would be produced by the benchmark physics models have different 
characteristics that are described below. 

The pruned jet mass and mass drop distributions in data, signal, and background simulations 
are shown in Fig. |2] The discriminating power of the pruned jet mass and mass drop for the 
different signals is evident. 

The product of signal B(W/ Z — > jets) and selection acceptance is shown in Fig.|3]for each model 
in the dijet invariant mass range up to 2 TeV (3 TeV) relevant for the singly(doubly)-tagged data 
analysis. The fraction of events that produce dijet events, which have survived the kinematic 
selection, is between 26% and 47%. This fraction includes also the branching fraction of W/Z 
decaying into objects which are reconstructed as jets, while it does not contain the W/ Z-tagging 
efficiency. The different behaviour of the acceptance for W and Grs at low dijet masses is due 
to different angular distributions generated by pythia 6 and HERWIG++. 

The W/Z-tagging efficiency for signal and background events is shown in Fig.|4] The signal ef- 
ficiency, determined from the simulation, is found to be between 20% and 45% (8% and 22%) for 
single (double) W/Z-tagged signals. The simulation modeling of the signal efficiency is cross- 
checked by measuring the W/ Z-tagging efficiency in semileptonic tt data, and by comparing 
it with the same efficiency obtained using the same procedure for tt simulation generated with 
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Figure 2: Pruned jet mass (left) and mass drop (right) in data, signal and background simu- 
lations. The signal simulation distributions are plotted as smooth curves connecting the his- 
togram entries (using the same binning as the data distribution). All simulation distributions 
have been scaled to match the number of data events. 
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Figure 3: Branching fraction into dijet final states B(W/Z — > jets) times angular acceptance 
(|//| < 2.5, \Arj\ < 1.3). The W/Z-tagging efficiencies are excluded from the acceptance. 

MadGraph 4.4.12 (sol and showered with pythia 6. The ratio of the two efficiencies results 
in a scale factor of 0.98 ± 0.03 which is then applied to the efficiencies for signals in the dijet 
data. The uncertainties on the scale factor are propagated into the systematic uncertainties on 
the overall signal efficiency. 

As described above, the production and decay of the Grs is modeled with HERWIG++. A differ- 
ence of up to 18% on the double-tag efficiency in the RS graviton WW/ZZ signal simulations 
between pythia 6 and HERWIG++ is observed. The difference can be attributed in equal parts 
to the different showering algorithms and the hadronization algorithms. The different under- 
lying event modeling has only a small impact of <1%. This discrepancy is accounted for as a 
systematic uncertainty on the double-tag efficiency. For the single-tag efficiency a 9% uncer- 
tainty is assumed. It is estimated from the double-tag efficiency uncertainty of 18% assuming 
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5 Signal characterization 
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Figure 4: Efficiency of requiring 1 W/Z-tag (left) and 2 W/Z-tags (right) in data, signal and 
background simulations for events passing the angular acceptance requirement (|//| < 2.5, 
\^r]\<l3). 
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Figure 5: The signal dijet invariant mass distributions for 1.5 TeV Grs — > WW, Grs — > ZZ, 
W WZ, q* qW, and q* qZ resonances, computed using anti-Zcj jets with R = 0.5. 
Crystal-Ball functions fit to the simulated distributions are shown. 

an identical difference in tagging efficiency for the two vector bosons in the case of 100% effi- 
ciency, representing an upper limit on this uncertainty. 

The effect of pileup on the W/ Z-tagging efficiency was also checked. Because of the rejection of 
charged particles not originating from the primary vertex and also the application of pruning, 
the pileup dependence is weak and the uncertainty of the modeling of the pileup distribution 
is less than 2%. 

The dijet mass dependence of the W/ Z-tagging efficiency for background events shown in 
Fig. |4] is adequately described by the simulation. Therefore, no additional systematic uncer- 
tainty is assigned on the dijet mass dependence of the modeling of the W/Z-tagging in simu- 
lation. 
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Figure [i] shows the signal shapes for Grs ZZ/WW, W' — > WZ, and q* — > qW/qZ, all of 
which correspond to a resonance mass of 1.5 TeV. The differences for the different models are 
due to the different tagging efficiencies for W and Z and the different modeling in pythia 6 
and HERWIG++. Resonance shapes were simulated at masses of 1, 1.5, 2 and 3 TeV and a linear 
interpolation was used to obtain the shapes at intermediate masses. 




Figure 6: The single (left) and double (right) W/Z-tagged mjj distributions (points) in data fit- 
ted with the QCD background parametrization (solid curve). For the double W/Z-tagged dis- 
tribution P3 = is assumed in Eq. ([l]|. Signal shape distributions for q* — > qW and Grs — ^ WW 
with arbitrary cross section are also shown. Bottom panes: the corresponding pull distributions 

/■ Data— Fit \ 



6 Background shape parametrization 

The shape of the QCD background in the dijet spectrum is modeled using a simple parametriza- 
tion which has been successfully deployed in previous searches in the dijet mass spectrum USTI . 
The background model is given by: 

d^ ^ Pojl-m/^sY^ 

dm (m/v/s)^2+^3ln(m/v^) ' ^ ^ 

where m denotes the dijet mass and a/ s the pp center of mass energy. Pq acts as a normalization 
parameter for the probability density function, and Pi, P2, P3 describe its shape. For the single 
W/Z-tagged analysis, all parameters are free to float in the fit. For the double W/Z-tagged 
analysis, P3 is not needed as suggested by a Fisher F-test and a simpler parametrization with 
P3 fixed to is used. 

Figure |6] shows the dijet mass spectra from single and double W/ Z-tagged fitted to Eq. ^ 
and the bottom panes show corresponding pull distributions, demonstrating the agreement 
between the background-only probability density function and the data. 

Since no sizeable deviation from the background-only hypothesis is seen, exclusion limits are 
set on the product of cross section, acceptance, and branching fraction for the five considered 
final states: qW, qZ, WW, WZ, and ZZ. 
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8 Limit setting procedure 



7 Systematic uncertainties 

The sources of systematic uncertainties are summarized as follows. The only background- 
related systematic uncertainty is the choice of background parametrization which is discussed 
in Section 8] The leading signal-related systematic uncertainties are the W/ Z-tagging efficiency 
(Section |5 1, jet energy scale (JES), jet energy resolution (JER), and luminosity measurement. 
Because the trigger and reconstruction efficiencies are larger than 99% in the relevant dijet mass 
range, the uncertainties associated with these efficiencies are negligible. 



In the jet pj- and )y-regions considered in this analysis, the JES has an uncertainty of 2-3% 
The Pt- and 7/-dependent uncertainty is propagated to the reconstructed dijet invariant mass, 
resulting in an uncertainty of 2.2%, which is approximately mass independent. The effect of the 
JES uncertainty on the calculation of the limits is estimated by varying the resonance dijet mass 
in the statistical analysis. The JER is known to a precision of 10% and its tails are in agreement 
between data and simulation Il44ll . The effect of the JER uncertainty on the calculation of the 
limits is estimated by varying the reconstructed resonance width in the statistical analysis. The 
luminosity has an uncertainty of 2.2% ||52| , which is also taken into account in the statistical 
analysis. 



8 Limit setting procedure 

For setting upper limits on the resonance production cross section a Bayesian formalism with 
uniform prior for the cross section is used, following the procedure used in Ref. [51] . The 
binned likelihood, L, can be written as: 

where 

m = aNi{S) + Ni{B), (3) 

n, is the observed number of events in the dijet mass bin, and N, (S) is the expected number 
of events from the signal in the i*^ dijet mass bin, a. scales the signal amplitude, and N;(B) is the 
expected number of events from background in the dijet mass bin. The background N, (B) 
is estimated as the background component of the best 5(4)-parameter fit of equation |3] to the 
singly(doubly)-tagged data points. The signal is not restricted to be positive for the background 
estimate fit although it is restricted in the Bayesian prior for the signal. A flat prior in a, which 
is the same as a flat prior in the resonance production cross section, is assumed. 

The dominant sources of systematic uncertainties (the jet energy scale, the jet energy resolution, 
the integrated luminosity, and the W/Z-tagging efficiency) are considered as nuisance param- 
eters associated to log-normal priors. The uncertainty on the background shape is taken into 
account with nuisance parameters associated to Gaussian priors representing variations of the 
fit parameters along the eigenvectors of their correlation matrix. The systematic uncertainties 
are accounted for using a fully Bayesian treatment and integrating the likelihood over nuisance 
parameters. 

The 95% confidence level (CL) upper limits Cs is calculated from the normalized posterior prob- 
ability density Ppost as follows: 

Ppost (cr)dc7 = 0.95 . (4) 
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This method of using the data first to constrain the background fit and second to extract the 
limit induces a bias in the coverage of the limits. The actual coverage is reduced to 93.6% 
(94.3%) at a WW (qW) signal mass of 1200 GeV (1800 GeV). 

9 Results 

Figure|7|shows the 95% CL cross section upper limits derived from the single and double W/Z- 
tagged event samples. The predicted cross sections as a function of resonance mass for the 
considered benchmark models are overlaid. A 95% CL lower limit is set on the mass of ex- 
cited quark resonances decaying into qW (qZ) at 2.38 TeV (2.15 TeV) whereas a limit of 2.43 TeV 
(2.07 TeV) is expected. These are the most stringent limits in the qW and qZ final states to 
date. The sensitivity of our measurement with the present dataset is not sufficient to extract 
meaningful mass limits on the Grs with fc/ Mpi = 0.1 nor the heavy SM-like W boson, but the 
cross section limits are the most stringent in the fully hadronic final state to date. Comparing 
to the cross section limits on Grs to ZZ in the corresponding semileptonic final states [20^-^22], 
this analysis sets a stronger cross section limit above a Grs mass of roughly 1.4 TeV. The cross 
section limits on final states with ZZ are stronger than those with WW because the efficiency 
for tagging a Z is larger than that for a W. On the contrary, the predicted cross section for WW 
is twice as large as for ZZ in the Grs model because of charge-conjugation. 

The dominant source of systematic uncertainty is the jet energy scale. Removing all systematic 
uncertainties on the single and double tagged searches would decrease the upper limit on the 
cross section by less than 20% at all values of resonance mass, which would translate into a 
change on resonances mass limits of roughly 3% . 

10 Summary 

A data sample corresponding to an integrated luminosity of 5.0 fb^^ collected in pp collisions at 
\/s = 7 TeV with the CMS detector was used to measure the W/Z-tagged dijet mass spectrum 
using the two leading jets within the pseudorapidity range < 2.5 and with pseudorapidity 
separation \Ati\ < 1.3. By suppressing the QCD background with the selection of candidates 
for vector bosons decaying to hadrons within jets, no evidence was found for new particle 
production in the W/Z-tagged dijet spectrum. A 95% CL lower limit is set on the mass of 
excited quark resonances decaying into qW (qZ) at 2.38 TeV (2.15 TeV) and upper limits on the 
cross section for resonances decaying to qW, qZ, WW, WZ, or ZZ final states. These are the 
most stringent limits in the qW and qZ final states to date. 
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Figure 7: Expected and observed limits for qW (top-left), qZ (top-right), WW (center-left), ZZ 
(center-right) and WZ (bottom) resonances. Here, B x A in the vertical axis label contains 
the branching fraction of Grs — > WW/ZZ — > 2 jets or q* — > qW/qZ — > 2 jets, as well as the 
acceptance for reconstructing the jets in < 2.5, \Arj\ < 1.3. The predicted cross sections as a 
function of resonance mass for the considered benchmark models are overlaid. 
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